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We use a cellular automata model to study the evolution 
of HIV infection and the onset of AIDS. The model takes into 
account the global features of the immune response to any 
pathogen, the fast mutation rate of the HIV and a fair amount 
of spatial localization. Our results reproduce quite well the 
three-phase pattern observed in T cell and virus counts of 
infected patients, namely, the primary response, the clinical 
latency period and the onset of AIDS. We have also found 
that the infected cells may organize themselves into special 
spatial structures since the primary infection, leading to a 
decrease on the concentration of uninfected cells. Our results 
suggest that these cell aggregations, which can be associated 
to syncytia, leads to AIDS. 

PACS: 02.70c,05.45-a,87.15a,87.18Hf 



The human immunodeficiency virus (HIV), which 
causes AIDS (the acquired immunodeficiency syndrome) 
has been the subject of most intense studies, that en- 
compass diverse fields of scientific research, ranging from 
the natural to the social sciences. Major progress has 
been achieved by medical and biological researchers in 
understanding the genetic code of the virus, the role of 
its high mutation rate, the virus-host interaction, the ap- 
parent failure of the immune system to control and elim- 
inate the virus, and the onset of AIDS. Nevertheless, the 
mechanisms by which HIV causes AIDS still remain un- 
explained. 

The immune response to any virus pathogen (virus, 
bacteria, etc.) is generated by a complex web of interac- 
tions, involving the cooperative and collective behavior 
of different types of white blood cells (monocytes, T- and 
B-cells). The time scale involved to develop a specific im- 
mune response may vary from days to weeks. The com- 
plex dynamic of HIV infection and the ensuing onset of 
AIDS involves a wide range of time scaleso. The primary 
infection exhibits the same characteristics as any other 
viral infection: a dramatic increase of the virus popu- 
lation during the first 2 — 6 weeks, followed by a sharp 
decline, due to the action of the immune system. For 
HIV, however, instead of being completely eliminated af- 
ter the primary infection, a low virus burden is detected 
for a long asymptomatic time: the clinical latency pe- 
riod. This period may vary from two to ten (or more) 
years, depending on the patient. Besides the low con- 
centration of HIV detected during this period, a gradual 
deterioration of the immune system is manifested by the 



reduction of CD4+T-cell populations in the peripheral 
blood. The third phase of the disease is achieved when 
the concentration of the T-cells is lower than a critical 
value (~ 30%), leading to the development of AIDS. As a 
consequence, the patient normally dies from opportunis- 
tic diseases, which would be usually controlled by the 
immune system. This common aattern observed during 
the course of the HIV infectionEl is depicted in Fig. 1, 
which shows the plasma viremia titer and the CD4+T 
cell counts in the Bcripheral blood as functions of time. 

Several theorie^ have been proposed to explain why 
and how the virus remain (albeit in low concentrations) 
in the organism after the primary immune response, and 
the causes of the decline of T-cell counts, leading to the 
onset of AIDS. So far, none of them has provided, a com- 
plete explanation for the entire process. The most com- 
mon assumption made to explain the permanence of the 
virus is that of fast mutation rate of the HIV, duCp, 
its rapid replication and its high degree of variability 
This assumption plays a central role in our theory as well. 

The challenge to understand the dynamics of the HIV 
infection has also attracted the attention of mathemati- 
cians and physicists that have been working on mathe- 
matical modelda to describe different aspects of the dy- 
namics of the host-parasite interaction. Although many 
of these models have contributed to the understanding of 
various aspects of the development of the disease, they 
fail to describe either the changes in the amount of virus 
detected in the body (blood, tissues and other body flu- 
ids) during the entire process or the two time scales ob- 
served in the course of the HIV infection, i.e., the primary 
response which runs from days to weeks and the clinical 
latency period that may vary from months to years. From 
the dynamical point of view, these different time scales 
may be related to two distinct kinds of interactions: one 
local and fast; and the other long-ranged and slow. 

Most of the mathematical models proposed so far use 
population equations, treating the infected patient as a 
homogenous entity. Therefore local interactions and spa- 
tial inhomogeneities, caused by localization of the initial 
of immune response in lymphoid organs, are not taken 
into account. We believe that this feature, which is a nat- 
ural ingredient of our model, are of central importance. 

Experimental evidenceaS support that the lymphoid 
tissue is a major reservoir of HIV infection in vivo. More- 
over a snapshot of the distribution of cells among the 
different compartments of the immune system will show 
only a small fraction (2 — 4%) of the cells circulating in 
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blood and lymph, while the majority is found in the lym- 
phoid organsQ. Paradoxically, the process of mobilization 
and activation of immune cells directed against the virus 
by the immune system, which occurs in the lymphoid 
micro-environment in this case_ provides a milieu that 
contributes to the virus spread.Q. 

Such spatial structure emerges naturally in discretg 
models, based on cellular automata that were shownQ 
to describe well cooperative and collective^ patterns in 
experimentally observed immune responses. Therefore 
we model the course of the HIV infection by a cellular 
automaton. Our model takes into account the main fea- 
tures of the immune response to any pathogen, the high 
mutation rate of the HIV and a fair amount of spatial 
localization that may occur in the lymphoid tissues. Our 
aim is to test whether the combination of these hypothe- 
ses can explain the three-phase dynamics observed on the 
course of the HIV infection (see Fig. 1). The results ob- 
tained by simulations of our model are shown in Fig. 2 
and, as far as we know, this is the first time that the 
complex dynamics of the HIV infection process has been 
so faithfully reproduced by a theoretical model. 
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FIG. 1. The common pattern exhibited by infected pa- 
tients, first presented by Pantaleo et altJ. The plasma viremia 
titer (black circles) and CD4^T cell counts (squares) versus 
time show a three-phase dynamics. 

We model the immune system cells in the lymphoid 
tissues by a two-dimensional cellular automaton. To each 
lattice site we associate a cell, that may be a target for 
the HIV (as for instance, a CD4-I-T cell or a monocyte). 
Each cell can be in one of four states: (a) healthy; (b) 
infected- Al, corresponding to an infected cell that is free 
to spread the infection; (c) infected-A2, an infected cell 
against which a specific immune response has already 
developed and, finally (d) dead, an infected cell that was 
eliminated by the immune response. 

The initial configuration is mainly composed of healthy 
cells, with a small fraction, pniv, of infected-Al cells, 
representing the initial contamination by the HIV. In one 
time step the entire lattice is updated in a synchronized 
parallel way, according to the rules described below. The 
updated state of a cell depends on the states of its eight 
nearest neighbors. 

Rule 1: Update of a healthy cell: (a): if it has at least 
one infected-Al neighbor, it becomes infected-Al. (b): if 



it has no infected-Al neighbor but does have at least R 
(2 < R < 8) infected- A2 neighbors, it becomes infected- 
Al (c): otherwise it stays healthy. 

Rule la mimics the spread of the HIV infection by con- 
tact, before the immune system had developed its specific 
response against the virus. Rule lb represents the fact 
that infected-A2 cells may, before dying, contaminate a 
healthy cell if their concentration is above some thresh- 
old. 

Rule 2: An infected-Al cell becomes infected-A2 after r 
time steps. 

An infected- A2 cell is one against which the immune re- 
sponse has developed and hence its ability to spread the 
infection is reduced. Here t represents the time required 
for the immune system to develop a specific response to 
kill an infected cell. Such a time delay is required for 
each infected cell since in our model we view each new 
infected cell as carrying a different lineage (strain) of the 
virus. This is the way we incorporate the fast mutation 
rate of the virus in our model. When a healthy cell is 
infected, the virus uses the cell's DNA in order to tran- 
scribe its RNA and replicate. During each transcription 
an error may occur, producing, on the average, one mu- 
tation per gefteration and hence a new strain of the virus 
is producedo'cl. 

Rule 3: Infected- A2 cells become dead cells. 

This rule simulates the depletion of the infected cells by 

the immune response. 

Rule 4- (o-) Dead cells can be replaced by healthy cells 
with probability Prepi in the next time step, (or remain 
dead with probability 1 — Prepi)- (b) Each new healthy 
cell introduced may be replaced by an infected-Al with 
probability Pinfec- 

Rule 4a describes the bone marrow replenishment of the 
depleted cells, mimicking the high ability of the immune 
system to recover from the immunosupression generated 
by infection. As a consequence, it will also mimic some 
diffusion of the cells in the tissue. Rule 4b simulates the 
entrance of new infected cells in the system that could 
come from other compartments of the immune system. 

This four-state automaton allows transitions from each 
of its four states into the next one in a cyclic manner. 
We performed simulations of the model on square lat- 
tices with periodic boundary conditions oi N = L"^ sites, 
with L ranging from 300 up to 1000. All the parameters 
adopted in the simulations were based on experimental 
data. The initial concentration of HIV, phiv = 0.05, was 
chosen based on the observation that the order of one in 
10^ or ll£. T cells harbor viral DNA during the primary 
infectionliy. We adopted Pmfec — 10~^ due to the fact 
that only one in 10^ to 10^ cells in the peripheral blood of 
infected individual express viral proteins. Since the im- 
mune system has a high ability to replenish the depleted 
cells, we used Prepi = 0.99, although smaller probabilities 
could also be considered, representing different individu- 
als. Since the time delay parameter (r) may vary from 
2 to 6 weeks, we chose t = 4. Each of our time steps 
corresponds to one week. 
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We present in Fig. 2 the evolution of the densities 
of healthy cells, infected cells (considering both Al and 
A2 types) and dead cells, obtained in simulations of our 
model . We show results averaged over one hundred sim- 
ulations (that used different initial configurations) and 
the corresponding standard deviations (shown as error 
bars). There is excellent qualitative agreement between 
our results for the density of healthy and infected cells 
and the time evolution of the number of CD4-I- T cells in 
the peripheral blood and the plasma viremia titer showed 
in Fig. 1. The model reproduces the two time scales 
observed in the dynamics of the HIV infection: a short 
one associated with the primary infection and a long one 
related to the clinical latency period and the onset of 
AIDS. The small error bars related to the first stage of 
the infection indicated that its dynamics is insensitive to 
the initial configuration. However, the large error bars 
during the latency period suggest the presence of some 
mechanism regulating this phase of infection. 
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FIG. 2. The results obtained from our simulations for a 

two-dimensional lattice with L = 700, Phiv ~ 0.05, R = 4, 

T = 4, Pin fee ~ 10~^, Prepi ~ 0.99. The evolutiou of the 

population densities exhibit the same three-phase dynamics 

observed for infected patients. We have adopted open squares 

for healthy cells, full circles for infected cells. 

The results presented in Fig. 2 are global properties, 
i.e. average quantities over the entire system. Our model 
allows also a closer look at local behavior, which in fact, 
may provide the clue for understanding also the global 
properties. From the analysis of the spatial configura- 
tions generated by the model in various individual sim- 
ulations, we noticed that the slow dynamics observed in 
the latency period is related to the emergence, after com- 
pletion of the primary response, of some special spatial 
structures of infected cells. These growing special struc- 
tures spread the infection in such way that they slowly 
commit more and more healthy cells, segregating and 
trapping uninfected cells. 

In Fig. 3 we present four "snapshots" of typical con- 
figurations obtained during one particular simulation. 
Starting form an initial configuration mostly composed 
of healthy cells (blue) with a random distribution of 
infected- Al cells (yellow), in subsequent time steps each 
individual infected-Al cell generates a pulse of infected 
cells, of width {t + 1), propagating in all directions. 
Whenever the average distance {h) between individual 



infected cells in the initial configuration is less than or 
equal (2r -l- 1), the independent pulses achieve a maxi- 
mum coverage of the lattice. For r = 4 the maximum 
coverage occurs after five weeks as shown in Fig. 3a. 
Note that the distribution of dead (red) cells corresponds 
to the initial configuration of infected cells. After that 
the concentration of infected cells decreases to a minimal 
value at 2(t + 1) time steps, establishing the end of the 
primary infection phase (10 weeks, in this case). 




FIG. 3. (Color) Four snapshots of parts of the lattice con- 
figuration for difi^erent time steps:(a)-(d) correspond to 5, 18, 
25 and 200 weeks, respectively. During the latency period 
we observe an organization of different types of cells in layers 
which interpenetrate each other. We have adopted the same 
parameters used in Figure 2. 

In the following time steps the presence of infected 
cells will be dictated by Pmfec, according to rule 4b. 
When new infected cells are introduced, they may gener- 
ate two different kinds of structures. The simplest case 
corresponds, as discussed above, to a wave of infected 
cells propagating in all directions, but in this case since 
Pin fee << Phiv the average distance {I2) between the 
new infected cells, introduced with probability Pinfee at 
random locations, is much larger than (2t + 1). Exam- 
ples of such structure are shown in the bottom right part 
of Figs. 3b and 3c. The spread of infection generated 
by these structures alone takes a longer time to reach 
nearly complete lattice coverage, and after that the in- 
fection vanishes. It is very rare, however, to find only 
these simple structures. A second type of structure also 
occurs, due to an interplay between rules lb, 4a and 4b. 
These special structures are generated by "sources" of 
infected cells. These sources appear, for instance, when 
a new infected-Al cell, introduced by rules 4a and 4b, 
is surrounded by at least four dead cells. At every ev- 
ery (r -|- 3) time steps they launch a propagating wave 
front of infected cells with width (r -I- 1). The period 
of the wave fronts corresponds to the period of time the 
infected-Al cell (source) remains in the same state (r) 
plus the three time steps to transition to types infected- 
A2, dead and healthy, in this sequence. Figs 3b and 
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3c (upper left) show the growing of such structures for 
two subsequent periods, corresponding respectively to 18 
and 25 weeks. As these structures grow, the number of 
infected cells increases and the concentration of healthy 
cells decreases. Moreover they segregate and trap un- 
infected cells between two consecutive waves of infected 
cells. These special pattern of cell aggregations explain 
the slow decrease of T cells in infected patients during the 
latency period. These growing structures may eventually 
cover the entire lattice, with final densities of different 
cells evolving towards a steady state with average con- 
centrations: {t + 1)/{t + 3) for all the infected cells, and 
1/(t -I- 3) for both the dead and healthy cells. Note that 
the average density for healthy cells at the steady state 
is always lower than the observed critical value(^ 30%) 
of CD4+ T cells counts related to the onset of AIDS in 
infected patients, i.e, the system's breakdown. Therefore 
the dynamics of real experimental data is related to the 
transient behavior of our model and not to its steady 
state. The long time scale, associated with this transient 
period that follows the primary infection, during which 
the model relaxes towards its steady state, corresponds 
to the clinical latency period. 

Analysis of the source distribution at any given 
time shows that the latency period depends on (^3), 
the average distance between sources and, conse- 
quently, on the probability of occurrence of sources (~ 
Prepi-Pinfec-PHiv)- Actually, since new sources can be 
released at any time step, the length of this transient 
time depends on the spatio-temporal average of the dis- 
tance between sources. 

In this work we have shown that our cellular automa- 
ton model reproduces quite well the three-stage dynamics 
observed in the course of the HIV infection, and the dif- 
ferent time scales observed in this process. The short 
time scale, characteristic of the primary infection, in- 
creases when T is increased or when phiv decreases. The 
long time scale, responsible for the clinical latency period 
and the onset of AIDS, is associated with the emergence 
of special structures, that increase slowly the number of 
infected cells and confine healthy cells. This special pat- 
tern formation depends on the value of the parameter 
R of rule lb, Prepi and Pmfec- Smaller values of Prepi 
enhance the presence of dead cells and favor the forma- 
tion of the special growing structures. We have also per- 
formed a careful investigation of the parameter space, 
and found that three-phase dynamics is reproducible for 
a wide range of the parameters. The complete study of 
the parameter space, the detailed discussion of the neces- 
sary conditions to generate the special spatial structures 
and the role they play spreading the infection, is under 
preparation and will be published elsewhere. 

The growing structures of infected cells may be asso- 
ciated with syncytia, the aggregation of cells observed 
experimentally, suggesting according to our results that 
they are responsible for the depletion of T cells leading to 
AIDS. These results actually corroborate some previous 
suggestions that syncytia are responsible for the perma- 



nence of the virus in the system, based on the analy- 
sis of the-similarities between HIV infection and other 
diseasetflO. 

Finally we emphasize that our cellular automaton 
model is, as far as we know, the first model which suc- 
ceeded to reproduce the general features of the complex 
dynamics of the HIV infection, as observed in infected 
patients. This woijk further substantiates claims made 
in previous studiesEI, that discrete models may be use- 
ful to describe emergent properties of complex biological 
systems, and to understand the mechanisms underlying 
its dynamical behavior. The reason for our success in 
describing the three-stage dynamics, whereas the other 
approaches fail, is that they do not take into account the 
spatial localization effects that play a major role on the 
course of the infection. A detailed study of the mecha- 
nisms underlying the dynamics described by our model 
may lead to new differential equation approaches, more 
suitable to describe the kind of dynamics observed in the 
course of HIV infection. 
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